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Abstract 

In the present work we investigate the behavior of a vortex in a long superconducting cylinder 
near to a columnar defect at the center. The derivations of the local magnetic field distribution 
and the Gibbs free energy will be carried out for a cylinder and a cavity of arbitrary sizes. From 
the general expressions, it considered two particular limits: one in which the radius of the cavity 
is very small but the radius of the superconducting cylinder is kept finite; and one in which the 
radius of the superconducting cylinder is taken very large (infinite) but the radius of the cavity 
is kept finite. In both cases the maximum number of vortices which are allowed in the cavity is 
determined. In addition, the surface barrier field for flux entrance into the cavity is calculated. 
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I. INTRODUCTION 

The analysis of vortex motion in the presence of defects is one of the most important 
topic in the superconductivity research field. The relevance of this theme is based on the 
fact that the behavior of type-II superconductors in the presence of an applied magnetic 
field can be described by the motion and pinning of vortices in the material. 

In the last decades a number of magnetic anomalies, such as the Wohlleben effect, fish 
tail anomaly and the jumps on the magnetic response of mesoscopic samples, were reported 
relating the occurrence of defects and the granular conformance of both bulk samples and 
thin films |l|, |2| . Simultaneously, several initiatives using the most different theoretical ap- 
proaches have been carried out to describe granular superconductors and their characteristic 
magnetic behavior p,|4|. Today it is well known that all these features and anomalies are 
straight linked to the vortex dynamics. 

Recently, the advances in the fabrication process related to nanotechnology made possible 
;he production of mesoscopic rings p and superconducting arrays with columnar defects 
6(. Experiments using these nanoscopic samples showed that the connectivity and vorticity 
are quite different than those for the macroscopic ones. The physics revealed by those 
experiments has contributed to a major understanding of type-II granular superconductors. 
This has motivated both experimental and theoretical physicists to investigate the magnetic 
properties of either small superconductors or bulk samples with nanosized defects. The 
properties of the vortex lattice, for instance, in this superconductors of confined geometries 
change radically with respect to the bulk superconductor ones. 

One of the open questions within the framework of mesoscopic superconducting materials 
is the validity and complete understanding of the simple rule proposed by Mkrtchyan and 
Shmidt pj. This rule indicates that the maximum number of vortices that enter a columnar 
defect is n s = r/2£, where r is the radius of the cylindrical cavity and £ is the superconducting 



coherence length. A 
cavity occupation 



though this rule has been questioned, it is a reasonable estimate of the 



In this paper the work of Refs. |2l llCj is extended for the case of a long superconducting 
wire, of arbitrary size, with a columnar defect. The applied magnetic field is parallel to 
the cylinder axis and the cylindrical cavity axis coincides with the direction of the applied 
field. First we determine the local magnetic field with appropriate boundary conditions 



and further, we calculate the Gibbs free energy of an ensemble of vortices near the cavity 
with n vortices inside it. On using the Gibbs free energy, the force acting on a vortex near 
the cavity and near to the external surface of the superconductor is calculated. From the 
expression of the force, we find a criterion for the maximum number of vortices allowed in 
the cavity. It is shown that if the size of the cavity is not too small, the saturation number 
differs substantially from the classical rule of Mkrtchyan and Shmidt 7]. The surface barrier 
field for flux penetration into the cavity is also determined. 

II. THE MAGNETIC FIELD 

The starting point of our study is the London equation for the local magnetic field of a 
very strong type-II superconductor for which the Ginzburg-Landau parameter k — j ^> 1. 
For our purposes, this equation can be more conveniently written in cylindrical coordinates 

as 
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where p t is the position of the i-vortex line, <J>o is the quantum flux, and N is the number 
of vortices outside the columnar defect. 

Consider a long superconducting wire under an applied field H parallel to the cylinder 
axis. We assume that there is a columnar defect at the center of the cylinder. Let a be the 
radius of the cylinder and r the radius of the cavity (columnar defect). This is a Dirichlet 
problem in which the boundary conditions are given by 

h(a,<p) = H, h(r,(p) = H , (2) 

where Hq is the magnetic field inside the columnar defect. Due to the symmetry of the 
geometry under investigation, the field in the cavity is uniform and parallel to the cylinder 
axis. The determination of Hq will be discussed later on. 

The solution of the London equation under these boundary conditions may be found by 
means of the Green's function method. We have for the local magnetic field 

N 

h(p, if) = $ Y^ G (P> ¥> P*> <Pi) + Mp) ' ( 3 ) 

i=\ 

where G(p,<p, p',ip') is the Green's function and hi(p) is the solution of the homogeneous 
London equation as if no vortex were present (Meissner state). By using Eqs. (J2J), we can 



easily determine the homogeneous solution, 

hl{p) = Ha ^ + H^\. (4) 

A (r,o) A (r,a) 

Here 

A m (x,y) = K m (x/X)I m (y/X) - K m (y / X) I m (x / X) , (5) 

where K m (x) and I m (x) are the modified Bessel functions of second kind. 
The Green's function can be expanded in a Fourier series as 

1 oo 

G(p,^pf, V f) = - £ e^-^g m (p,p')- (6) 

m=— oo 

Although easy, it is very cumbersome to determine the Fourier coefficients g m {p, p')- Here, 
we will omit this calculation. The Dirichlet boundary conditions (J2J) impose that 

9m(r, p') = g m (a, p') = , (7) 

at the surfaces, and 

9m((p'T,P') = 9m((p') + ,p) , 

dg m \ _ fdg m \ 1 , . 

The first condition represents the continuity of g m function at p — p'. The second condition 
represents the discontinuity of the derivative of g m at p = p', which comes from the delta 
function in Eq. (IT1). By using these boundary conditions, we find 
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where A m (x, y) is given by Eq. (JSj). Notice that this function is antisymmetric by inter- 
changing x and y, that is A m (x,y) = —A m (y,x). 

We still have to find the constant H . For that, we use the supercurrent expression 

Vxh = i(l v "- A )- < 10 > 

where A is the vector potential, and 9 is the phase of the order parameter. By integrating 
this equation in the circle of radius r containing the cavity and using the fact that the phase 
is not single valued, we obtain 
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where n is the number of vortices which were captured by the columnar defect. It should 
be clear the difference between N and n. 

Only the m = term will survive in the integration of Eq. (jllj) . A length algebra leads 
us to the following expression for the magnetic field in the cavity 
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where 



S m (x,y) = K m+1 (x/X)I m (y/\) + K m (y/X)I m+1 (x/X) 

Notice that in the limit of large cylinder, a — > oo and r finite, A (r, a) 
and 6 (r,a) — *■ Ki(r/X)I (a/X). Then, we have 
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which is the result found recently in Ref. 
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III. THE GIBBS FREE ENERGY 



In order to proceed we need to evaluate the Gibbs free energy (per unit length) which is 

given by 

a r> u 

(15) 
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where F is the London free energy (the Helmholtz free energy in the thermodynamic context) 
and B is the average induction; here A is the area of the cylinder cross section. In the London 
approximation, the free energy is given by 
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By using integration by parts, we can write 
r $ o ^ . , v X 2 aH f 27r fdh\ 
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Next, we compute the derivatives and insert into this last equation. In addition, we use 
Eqs. © and (JU). We find 
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The average induction may be obtained by integrating the London equation (JIJ on the 
cross section. We have 
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Finally, upon substituting Eqs. (fT£J) and (fTHjl into Eq. (fTKJ) . we obtain for the Gibbs free 
energy 

8tt 



<J)^ " (J> E7 _, 

G = ^ ^ C(Pi, ¥>i, Pj, <Pj) + J-A? r J2[A (pi } a)+A (r, Pi )] 



2rr2 



+ 



A 2 # 



+ 



47rA (r,a) 
Arif 5 (r,a) ( H { 



i=X 



'0 



Aa# 2 <Wa,r) N$nH 



2 A (r,a) ' 2 A (r,a) V 2 7 4 A (r,a) !;r 

In the Appendix we show that the Green's function can also be written as 
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A m (r,a) 

where p> and p< are defined in Eqs. (jA2|) and ()A3|) . Notice that for p = p', London theory 
requires a regularization procedure. We use a sharp cutoff in which \p — p'\ is replaced by £. 
In some special cases in which we apply the above results, Eq. (}2*T]) will be shown much 
more useful than Eqs. (|BJ) and (JHJ). All the results above are general, so it should be applied 
to a superconducting cylinder with a columnar defect, both with arbitrary sizes. 



IV. FORCE 



We will be most interested in the particular case N = 1. From Eq. EH we have 
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Now the Green's function is given by 
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and the magnetic field in the cavity is 
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The force acting on a single vortex can be found by taking the derivative of the Gibbs 
free energy. At the surface of the columnar defect, the force (per unit length) is then given 
by 
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where the use of the identities Ki(x)Io(x) + K (x)Ii(x) = 1/x, K' (x) = —Ki(x), and 
Iq(x) = h(x) have been made. 

Following a similar procedure, we can find an expression for the force acting on a single 
vortex at the external surface of the superconductor. We have 
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The derivative of the Green's function can be easily evaluated on using Eq. (J2*3*j) . Very 
close to the cavity, we can write 
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Notice that at p! = r + £ the Green's function vanishes as required by the boundary condi- 
tions. Consequently, on using a sharp cutoff we have 
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where it has been used the approximation Ki(£/X) = A/£ which is valid for a strong type-II 
superconductor. 

Similarly, we can find for the derivative of the Green's function at the surface of the 
superconductor, 
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Notice that these terms of both forces do not depend on the geometrical factors, but only 
on the fundamental lengths. 

From now on we will take two special cases. One case will consist of a columnar defect 
very small compared to the penetration length A but still keeping a/A finite. In the other 
case we will take r/X finite and a very large cylinder, a/A — ► oo. 



A. Small Cavity, Finite Superconducting Cylinder 

We have experienced that to find a closed expression for the force, both at the cavity 
and the external superconductor surface, is very cumbersome, even in the case r/X <C .. 
Nevertheless, by using symbolic mathematics on computers this task becomes feasible |ll|. 
From Eqs. ()28|) and (|30|). up to second order in r/X, we find 



$ 2 1 



^vr 2 A 2 r [2£ 



(n + l) + 



~fr 



Io(a/X) 
where C = 0.577215 ... is the Euler's constant. 
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Upon neglecting terms of second order, then the force will be negative or vanish if n + 1 < 
r/2£. Therefore, the maximum number of vortices permitted into the cavity is given by 



n s =n + 1 = — . (31) 

This result was firstly found by Mkrtchyan and Shmidt [7J. However, their derivation 
was valid only for infinite superconductors. It is remarkable that for finite superconductors, 
this saturation number does not depend on the radius of the superconducting cylinder but 
only of the size of the cavity. So, even in the case of a mesoscopic superconductor, in which 
r<ttCA, the proximity of the external and internal surfaces has no influence on how the 
cavity will be occupied. 

Let us now turn our attention to the surface barrier field sufficient to saturate the cavity 
with vortex entrance. From Eqs. (|2fi|) and ()29j) . by keeping only terms of second order or 
less in a/ A, we find 
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If —f a > the vortex enters the sample. Then, the minimum applied magnetic field for 



entering n + 1 vortices in the sample is 
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where H%° = $o/4vrA^ is the surface barrier field of a bulk superconductor, and we have 
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used Eq. (jUTJ). Notice that for a/A > 1, I m (a/X) = e a/x /^/2ira/\. Thus, for a large 
superconductor, Eq. (J3~3~j) reduces to H s = H%° which is a well known result for a bulk 
superconductor [12j. On the other hand, if a/ A <^ 1, then J (a/A) = 1, I\(a/\) = a/2 A, so 
that 

a \2i 

The term inside the square brackets is of order or larger then the unity, so this field 
is much larger than the surface barrier field for a bulk superconductor. As expected, a 
mesoscopic superconductor would required a much larger applied field for flux penetration 
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B. Finite Cavity, Large Superconducting Cylinder 



Let us move on to the discussion of another special case in which the size of the cavity is 
finite and the size of the superconducting cylinder is infinite. Within these limits, the free 
energy of Eq. (|22j) becomes 
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where now 
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Thus, the force acting on a single vortex at the surface of the cavity is given by 
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The (n + l)-th vortex will be captured by the cavity if this number makes the force zero. 
Therefore, the saturation number is 
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Notice that for the case of a small cavity this last expression yields the Mkrtchyan-Shmidt 
rule. In Fig. ^we plotted Eq. (|38|) as a function of r/2£ for several values of k. We can see 
that, as the size of the cavity increases, the saturation number differs significantly from the 
Mkrtchyan-Shmidt classical rule. Thus, more vortices should be allowed in the columnar 
defect if the radius of cavity increases, however the larger k the lower the increasing is. 



V. SUMMARY AND CONCLUSIONS 

In summary, we have derived the local field and the free energy of an ensemble of vortices 
around a columnar defect in a superconducting wire for both the cavity and the supercon- 
ductor of arbitrary sizes. We also evaluated the force near to the cavity and external surfaces 
of the superconductor. It has also been found that for a large superconductor, the cavity 
saturates at a larger number of vortices as its size increases. However, not linearly in r/2£ 
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as predicted in pj. In fact, in [8|, |9| it was presented some numerical simulations which 
states that the population of the columnar defect grows not linearly with r/2£, although in 
a different manner as found in the present work. Perhaps because in the present work it 
has been used the London theory which is strongly dependent on the cutoff model for the 
vortex core, whereas in those references the Ginzburg-Landau theory was employed. 

20 




FIG. 1: The straight line is n s = r/2£. The other curves are given by Eq. IfMHJI : from the top to 
the bottom k = 20, 30, 40, 50. 



APPENDIX A: GREEN'S FUNCTION 



The expression found for the Fourier coefficients Q may be written in a more compact 

form as 

A m (a,p>)A m (p<,r) 
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and 
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Introducing the expressions of A m of Eq. (JHJ), we find 

9m(p,p') = X2A , r[K m (a/X)I m (r/X)I m (p > /X)K m (p < /\) 
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+K m {r/X)I m {a/X)K m {p > /X)I m {p < /X)} . (A4) 

Inside the square brackets we can sum and subtract the following term 
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Substituting this equation in Eq. © and using the identity 
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On using the definition (J3J) we arrive at Eq. (J2TJ). 
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